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ABSTRACT 


Radiation of energies from 10 keV to greater than 10 MeV has been observed 
during solar flares, and is interpreted to be due to bremsstrahlung by relativistic 
electrons. A complete treatment of this problem requires solution of the kinetic 
equation for relativistic electrons and inclusion of synchrotron energy losses. Using the 
elctron distributions obtained from numerical solutions of this equation the bremss- 
trahlung spectra in the impulsive X-ray and 7 -ray regimes are calculated, and the 
variation of these spectral indices and directivities with energy and observation angle 
are described. The dependences of these characteristics of the radiation of changes in the 
solar atmospheric model, including the convergence of the magnetic field, the injected 
electron spectral index, and most importantly, in the anisotropy of the injected electrons 
and the of convergence of the magnetic field are also described. The model results are 
compared with stereoscopic observations of individual flares and the constraints that 
this data sets on the models are discussed. Subject headings: Sun:flares-Sun:X-rays-X- 
r ay s : spectr a-r adiat ion mechnanisms 
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I. Introduction 


It is generally accepted that the hard X-rays and continuum 7-rays from impulsive 
solar flares are due to bremsstrahlung radiation from non-thermal electrons. The 
characteristics of the radiation such as the spectrum, directivity, polarization and spatial 
structure depend on the properties of the accelerated electrons and the geometry of 
the magnetic field. Consequently observed radiation can give information about the 
distribution of the accelerated electrons and the magnetic field structure during the 
impulsive phase. For this purpose we have developed programs to calculate the radiation 
expected from models with different distributions of accelerated electrons and magnetic 
field structures. The comparisonsof the model results with observations constrain the 
models. 

Using a steady state Fokker-Planck equation described by McTieman and 
Petrosian (1989), hereafter referred to as Paper I, we have calculated for many models 
the electron flux F(F,/i,s), where E is the kinetic energy, /i is the pitch angle cosine, 
and s is the distance from the injection point along the field lines (cf. Paper I). In this 
paper we shall describe the calculation of the spectral characteristics of the bremsstrah- 
lung radiation from the non-thermal electrons. This work is an extension of earlier work 
by Leach and Petrosian (1983) on X-ray radiation to emission of 7-rays by relativistic 
electrons. 

Observations of impulsive phase radiation in regard to its directionality have been 
difficult and the results somewhat controversial. There exist two kinds of observations. 
The stereoscopic observations carried out by Kane and his collaborators using the 
ISEE-3 (or ICE) and PVO satellites provide the most direct method of determining 
the isotropy of impulsive hard X-rays and 7-rays from 100 keV to 1 MeV, (Kane 
et al. 1988). However, at present there have been very few observations of strong 
flares from two largely separated directions. The second kind of observations relies 
on statistical variation of radiation characteristics such as rate of occurrence, burst 
intensity and spectral characteristics with flare position on the solar disk (i.e., observed 
disk to limb variations). The first such study, using observations from the OSO-7 


2 


satellite, was described by Datlowe et al. (1974) and analyzed by Petrosian (1975) and 
Datlowe et al. (1977). These observations of low energy hard X-rays show little evidence 
of anisotropy in agreement with the theoretical models of Leach and Petrosian (1983). 
The more recent observations from SMM also show similar results in the X-ray 
regime. This picture changes gradually with increasing energy; the anisotropy of the 
radiation appears to increase leading to the dramatic concentration of flares with 
> 10 MeV emission near the solax limb (Riger et al. 1984). The statistical behavior 
of the anisotropy in the 7 -ray regime obtained from GRS on SMM is described by 
Vestrand et al. (1987) and Bid (1988). This gradual change is evident in most of the 
model results and, as pointed out by Petrosian (1985), can be understood in terms of 
the decrease with energy of collisional pitch angle diffusion. 

In this Paper we concentrate on spectral variation with the viewing angle of 
non-thermal thick target models. The calculation of the bremsstrahlung radiation is 
discussed in Section II along with the parameters of the assumed model atmosphere and 
field geometry of the flare region. In Section III we describe the variation of the emitted 
radiation with observation angle, position within the flare loop and photon energy for 
models with different accelerated electron energy and pitch angle distributions, and 
different degrees of convergence of the magnetic field. The results are compared with 
the stereoscopic observations of flares presented by Kane et al. (1988) in Section III. 

In Section IV we summarize our results. In a third paper of this series we compare the 
model results with the observed disk to limb variations. 

II. Calculation of Bremsstrahlung Emission 

In order to get the radiation from an electron distribution we must integrate the 
product of the electron flux with the bremsstrahlung cross section. We follow the 
procedure described in Leach (1984) to evaluate the Stokes parameters, which are 
extensively discussed in the literature. [See e.g., Haug (1972) or Bekefi (1966).] The 
Stokes parameters are then used to calculate the spectrum and polarization of the 
radiation. 

The radiation is linearly polarized due to the reflective symmetry around the 
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plane of the photon and incident electron momenta. In this paper we will describe the 
spectrum and directivity but not the polarization of the radiation. The photon flux as 
a function of the photon energy k, and polar and azimuthal angles 0 and <f> is obtained 
from 


= n(s) I™ dE [*dQ I' d*<7(k,E,Q)F(E,v,s) , (1) 

JO Jo J —tt 

where n(s) is the density of nuclei and ions, and 0 and $ are the polar and azimuthal 
angles of the electron momentum with respect to the photon direction k. As described 
in Paper I electrons are tied to the field lines so that the electron flux may be expressed 
in terms of the electron energy E, pitch angle cosine ft and distance along the field line 
s. The cross section a is the stun of the cross sections for parallel and perpendicular 
polaraizations. It will not be reproduced here, but may be found in Koch and 
Motz (1959), equation (2BN). 

The electron pitch angle cosine ft is a complex function of the angles 0, 6 and <f>. 

For comparison with observations it is convenient to express the photon flux relative to a 
coordinate system fixed at the surface of the sun with the z sods in the radial direction. 
The polar angle 0 of the photon then corresponds to the heliocentric angle of the flare; 

0 = 0 for a flare at the center of the disk, and 0 = 90° for a flare on the limb. The pitch 
angle cosine is given by 


ft as cos$ e cos0 B + sin 0 e cos <£ e sin , (2) 

where 0 B is the angle of the magnetic field with respect to the 2 - axis, 


cos 0 e — cos 0 cos 0 — sin 0 sin 0 cos $ 


and 


<f> e =s <f> + tan -1 


sin 0 sin $ 

cos 0 sin 0 cos $ + sin0cos0j 


( 3 ) 


Given the variation along the magnetic field line of the density n(s) and the 
variation of the angle of the magnetic field with respect to the z axis, 0 B (s), then the 
flux in any direction may be found as a function of distance s from the above equations. 
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For the field variation we assume a magnetic loop consisting of a semicircular coronal 
portion with vertical footpoints below the transition region. 

Since the available observations of directivity are spatially unresolved we need the 
integral over depth of the photon flux, given by 

J(k,6, (j>) = (47rr^ In A) -1 / dE f dr c Td© f" d*<rF(E,fi,T c ) , (4) 

Jo Jo Jo J-it 

where r c is the dimensionless column depth defined in Paper I; dr c = 47rr^n(s)ln Ads, r 0 
is the classical electron radius, and In A « 20 for a fully ionized ambient plasma. 

Most of the radiation, especially at higher energies, comes from the footpoints below 
the transition region where the magnetic field is vertical. Hence most of the radiation 
will be azimuthally symmetric and depend only on the polar angle. The small amount of 
low energy radiation from the coronal part of the loop could be azimuthally asymmetric 
and could affect the center to limb variation of the radiation, depending on the position 
and orientation of the loop with respect to the solar equator. For example, for a loop 
straddling the equator, and the angle of the projection of the loop on the disk is 7r/2, 
there is no azimuthal asymmetry throughout the loop, while if the angle is less than 7t/2, 
the emission from the coronal portion will be asymmetric and the overall center to limb 
variation will be slightly different than that of the vertical part. In what follows we shall 
ignore this small dependence and use the azimuthally averaged photon flux 

J ( M ) = ( 2')- 1 

Finally we define the average flux 

J(k,9)sin0d6 (6) 

which will be useful especially for low energies where the bremsstrahlung emission is 
nearly isotropic. 

The evaluation of the electron flux F(E, fi,T c ) and the integrations in equations (1) 
said (4) require, in addition to $ B (a) and n(a), the knowledge of the variation of the 
strength of the magnetic field (actually the logarithmic derivative dlnB/ds) and plasma 



[ 2w J(k,e,t)d<f> . 
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density n(s) along the field lines. These sure obtained from the solar atmospheric models 
which will be the same as those discussed in Paper I. In Figure 1 we show the variation 
with hieght of the density n, temperature T and column depth N from the top of the 
loop. Note that the atmospheric model prescribes variation with hieght, not 3. Below 
the transition region, because of our assumption of vertical field lines, the two variations 
are identical. This is true for density and temperature in the corona as well because 
these quantities are constant there. The variation with s of column depth, however, 
depends on the radius of the loop or the length of the loop to the transition region s tr . 

We assume a pure hydrogen atmosphere; the addition of heavier elements 
renormalizes the results but has no effect on the angular or spectral variation of the 
radiation. The level of ionization of hydrogen is computed using the Saha equation and 
the temperature profile given in Figure 1. In the corona the plasma is fully ionized. The 
ionization level decreases with temperature, and the plasma is mostly neutral below the 
transition region. As shown in Paper I, the level of ionization can have some effect on 
the angular distribution of the electron flux so it might affect the angular variation of 
the radiation. 

In Figure 1 we show two possible variations for the magnetic field. In one it is 
assumed, as in Leach and Petrosian (1983), that Sy = dlnB/ds is constant, so 
that most of the convergence of the field occurs in the corona. Since most of the 
high energy radiation comes from the footpoints, from the electrons that escape the 
corona, we call this a “precipitation” model. In the second kind of model we assume 
dlnB/dlnn = constant, so that the field is constant in the corona but increases as 
some power of the density (or pressure, since T is nearly constant) below the transition 
region; B oc n*. This kind of convergence, which we call the “trap” model, may have a 
more significant effect and is similar to the models of MacKinnon and Brown (1988) and 
Miller and Ramaty (1989). 

We note, however, that a complete knowledge of the variation of B(s) and n(s) is 
not always needed. As can be seen from equation (4), for analysis of spatially unresolved 
observation the knowledge of n(s) is not necessary. Similarly, as shown in Paper I, for a 
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flare the synchrotron losses are important only for very high energies and if the magnetic 
field in the corona exceeds 10 4 G. If we ignore synchrotron losses we need only to specify 
the value of din B/dr e for a complete description of the problem and not B(s) and 
n(s ) separately. For precipitation models most of the convergence is in the corona and 
whether it influences the electron distribution depends on whether b T = B(s tr )/B 0 is 
greater than one. We specify this function which is directly related to the parameter 
dlnf?/dr e given the coronal density n 0 = 10 lo cm” 3 and the depth of the transition 
s tr = 2.4 x 10 4 km. For the “trap” models, the most important paramter is t but 
in practice, since the numerical integration can only be extended to a finite column 
depth iV m4X , we specify the ratio B(N mmx )/B 0 from which we evaluate e. We choose 
N ma . x = 6 x 10 25 cm” 2 which corresponds to r e « 650 at which point the vast majority of 
electrons with energies < 100 MeV have disappeared. 

As in Paper I we assume injected fluxes with power law energy spectra and isotropic 
or gaussian pitch angle distributions, with maxima at pitch angle c*i (<*1 = 0 for “beam” 
models or a\ = ir/2 for “pancake” models) and dispersion a 2 : 

F(£, /x, 0) = 2 tffr*a; 2 « _(a_ai)J/o ’ • (7) 

where small values of or 2 imply highly beamed or flattened distributions and a 2 — » 00 
implies an isotropic distribution. The total number of electrons with energy greater than 
E 0 is then KEl~ s /(6 - 1). 

III. Characteristics of Radiation 

The characteristics in which we are interested here are the directivity (variation 
with 9 and <f>) and the spectrum (variation with k) of the bremsstrahlung radiation. 

For comparison with spatially resolved observations we need in addition these 
chararcteristics as functions of depth s along the field lines, namely J(k,9, 4>,s). There 
are, however, very few reliable spatially resolved observations, especially at high 
energies. Spatially resolved observations by HXIS or Hinotori were primarily at low 
energies and sometimes suffered from contamination from the thermal component. 

The only information on the spatial structure at higher energies has come from 
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limb occulted flares (Kane 1982) which were analyzed by Brown et al. (1981) and 
Leach and Petrosian (1983). Consequently, we shall leave a complete description of 
depth dependence of emission for all models to the future if and when data becomes 
available from Max ’91 missions. For now we concentrate on the spatially integrated 
characteristics and their comparison with existing observations. However, in order to 
clarify some aspects of the spatially integrated results we present the characteristics of 
emission as functions of depth for one model. 

A. Dependences on Depth 

For this purpose we select the model, which will be called the standard model, with 
an isotropic injected pitch angle distribution (a* = oo) and a uniform magnetic field 
(dlnB/ds = 0). We shall be concerned with field strengths and energies such that 
the synchrtron losses are negligible and the problem can be formulated in terms of the 
column depth r c = N/N 0 , where N 0 — 5 x 10 22 . As shown by Leach (1984) (see also 
Canfield et al. 1985 or Brown and McClymont 1975) for non-relativistic energies the 
dependence of the bremsstrahlung intensity on depth scales very simply with energy. He 
found that the fractional emission integrated over all angles, I(k,r e ) = J(k,r e )/ J(k), 
is proportional to (Jb 2 + t c )~ 6 ! 2 . The logarithmic plots of I versus r c are very similar 
to the plots of the electron flux integrated over pitch angle F M (.E,r c ) versus r c shown 
in Figure 5 of Paper I. In Figure 2a we show the variation of r c I with r e . We have 
multiplied I by r c in order to show from where, on this logarithmic plot, most of the 
emission originates. As evident from the figure this happens when the curves reach their 
maxima at r c « rj(k) = k 2 /(k + 1). The simple scaling found by Leach (1984) gradually 
breaks down as energy increases. 

The actual spatial dependence of the emission depends on the variation of the with 
depth of the density (for the standard model) and the magnetic field (for the general 
models). To demonstrate this, in Figure 2b, we show the variation of the fractional 
emission I(k,s) = J(k,s)/J(k) with depth s. Since drjds <x n(s) these curves behave 
like the product of /(Jfc,r c ) and the ambient density. In the corona, where density is 
constant and r c <C k 2 /(k 4- 1), I(k,s) is constant. In and below the transition region, 
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the density increases rapidly, so I(Jc,s) increases as the density, but eventually beyond 
r c = k 2 /(k + 1), I(k,s) begins to decrease because I(Jc,t c ) decreases faster than the 
increase in density. Thus there is a maximum for fractional emission I(k, s) at the depth 
corresponding to Tc ss k 2 /(k + 1), which increases with energy. For N tT — 2.4 x 10 19 cm -2 
as assumed here, for energies greater than ~ 12 keV, the maximum in emission 
occurs below the transition region, so that the projected image will show two emission 
regions around the footpoints of the loop. Only for higher column depths N tt and/or 
lower energies will there be substantial emission from the coronal part of the loop, in 
which case the geometry and orientation of the loop are important in determining the 
projected image of the emission. For low coronal density or high energies (gamma rays 
> 300 keV) these complex geometric factors are absent and the bremsstrahlung emission 
has azimuthal symmetry so that nnnnthe spatially resolved or integrated emissions 
depend only on the polar angle 9. We define the total and depth dependent directivities 
as 

A (M) = J(k,0)/J(k,9 0°) , A (M,r c ) = J(M,r c )/J(*, 90 V c ) , (7) 

Note that 6 = 90® corresponds to a flare at the limb while 0 = 0° corresponds to a flare 
at the center of the disk so that the directivity is normalized at the limb. The directivity 
of the standard model at different depths for 1 MeV emission is shown in Figure 3. 

First note that even for isotropically injected electrons most of the radiation is directed 
towards the photosphere (90 < 6 < 180°), and only & small fraction is directed into the 
observable hemisphere (0° < 9 < 90°), especially at larger depths. In general, for a given 
model, the qualitative behavior of the directivity is the same as that of the electron 
pitch angle distribution, but the anisotropy in the radiation distribution is smaller than 
that for the electron distribution, since the cross-section has a finite width. 

Figure 4 gives the normalized photon spectrum J(k, versus energy k for 

limb flares (9 = 90°) at different depths. The curves labeled corona, Photosphere 
and 90 represent the spectra for emission from the corona, photosphere and whole loop, 
respectively. In the corona the spectra are steeper because the higher energy electrons 
emit mostly at lower depths. The spectra become flatter with increasing depth, and at 
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the photosphere the spectral indicies are less than <5—1, which is the appropriate value 
of the spectral index for the emission from the whole loop. 

. B. Characteristics of Unresolved Emission 

The spatially unresolved observations can be dealt with by consideration of the 
variation of the emission from the whole loop with photon energy k and observation 
angle 9 which is given by the flux J(k,9 ) and the directivity A (k,9). Te general 
behavior of these quantities is demonstrated by the dashed lines marked 0 and 60 and 
the solid line marked 90 in Figure 4. The spectra become flatter going from 9 = 0° to 
9 = 90° (center to limb) but are all steeper than the curve for the total emission (the 
dashed line labeled Total), which includes emission for angles greater than 90°, directed 
into the sun. This is a consequence of the change in anisotropy with energy depicted 
in Figure 5 which shows A(fc, 9) versus 9 for energies 16 keV, 300 keV, 1 MeV, and 
10.6 MeV. As expected the range of the directivity increases with energy. Most of the 
variation in the directivity occurs between 60° and 90°. The statistical evidence for such 
changes of spectral index with angle has been observed (Vestrand et al. 1987). Since the 
magnitude of the change is model dependent, such observations should enable us to put 
constraints on models. This will be the subject matter of our third paper of this series. 

The spectra can not be described well single power laws since the curves steepen 
or flatten as the energy increases. We describe this behavior by fitting the photon flux 
to power laws in three regimes and obtain the spectral indices y x for 16 to 210 keV, j m 
for 300 keV to 1 MeV, and 7 7 for 2.6 to 21.2 MeV. The breaks in the spectra, i.e., the 
differences between the spectral indices in different energy ranges, are also dependent on 
angle, and the way they change with angle is model dependent. We will discuss these 
dependences in the following section. 

Note that the total spectrum integrated over all angles has indices y x « 6 — 1, 

7 m « <5 — 1 and 7 7 < 6 — 1; the total spectrum flattens slightly for extreme relativistic 
enegies. This is generally true for all of the models, and is due to the fact that collisions 
are the dominant loss process. We can find this result easily using the non-relativistic 
and relativistic limits for the bremsstrahlung cross section integrated over all angles 
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[equations 3BN(a) and (b) from Koch and Motz (1959)], and by assuming that the 
electron flux is completely isotropic with the energy dependence given by equation (28) 


of Paper I; i.e., 


F tot (E) — oc 


£-(S-2) 
£-( 6 - 1 ) 


for E <C m e c 2 
for E m e c 2 


(8) 


We may then integrate the product of flux and cross section over energy, and we find 


'TOT 



k -(6-l) 

k -(S-i) In (fc/m e c 2 ) 


for k < m e c 2 
for k m e c 2 


(9) 


The relativistic regime has a slightly flatter spectrum due to the logarithmic dependence 
on energy. We now consider the effects of variation of the model parameters. 


1. Plasma Density Profile and Ionization Level 

The ratio of bremsstrahlung emission rate to collisional energy loss rate is fairly 
independent of the density and varies slightly with the ionization state of the plasma. 
Consequently, when collisions are dominant (synchrotron and other losses are negligible) 
the emission form unresolved sources is insensitive to the form of the density or 
ionization profile. This is demonstrated in Figure 5 where the directivity at 1 MeV 
of the standard model is compared with two models. The dashed line represents the 
directivity for the model with constant density and full ionization and the o’s represent 
the model with a density profile as shown in Figure 1 which increases from 10 lo cm -3 in 
the corona to 10 18 cm -3 at s = 4300 km below the transition region, but is fully ionized 
throughout. There is little difference among these curves, especially for 0° < 9 < 90° 
which is the range of observable radiation. 

2. Electron Spectral Index 

The angular variation is also fairly insensitive to the value of the injected spectral 
index 6. In general the anisotropy of the radiation increases with decreasing values of 
the spectral index of the electrons, but this increase is less than 10% for relavant values 
of 6 (approximately from 2 to 7). The photon spectral index, of course, varies with the 
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electron spectral index, but the variation with 6 of the difference between these two 
indices, 

T, = 7. - <* (i = x,m, 7) (10) 

is in general small. 

For the remainder of this paper we shall ignore these differences and, assuming 
6 = 3.75, discuss the effects of the electron pitch angle distribution parameters aj , 
and the convergence of the magnetic field (the parameter d\nB /ds or the mirror ratio 
b T ). For each set of parameters we show the dirctivity normalized to unity at 9 = 90° 
(limb) for observable angles (0° < 9 < 90°) and three representative energies, 16 keV, 

1 MeV and 10.6 MeV. Similarly we show the variation of the three spectral indices 
r„ r m and r 7 , as defined in equation (10). As we shall see, the general feature of the 
curves is a hardening (flattening) of the spectra with increasing energy and observation 
angle. 


3. Electron Pitch Angle Distribution 

The effects of the injected pitch angle distribution are shown in Figures 6a and 6b. 
We give results for three models spanning the possible range for the distributions, one 
highly beamed model (o^ = 0.04, ou = 0, dashed lines), one isotropic (a^ — ► 00, solid 
lines) and one highly flattened pancake model (a* = 0.04, oti — tt/ 2, lines). As 
in the standard model, the total anisotropy A(fc,0°), defined as the ratio of emission at 
0° to emission at 90°, and the total variation in spectral index, AT = r(0°) — r(90°) 
increase with energy in all cases. These increases are smallest for the beam injection 
(dashed lines), and become larger for the isotropic and pancake models model (solid 
lines and -+ lines, respectively). In each case the gradients are small between 0° and 60° 
and most of the variation takes place near the limb between 60° and 90°. 

As pointed out earlier (Petrosian 1985), this behavior is a consequence of the 
shape of the cross section, which in general has a maximum at a photon-electron angle 
0 = 0 of width A0 » m t <? /E, but is fairly flat at large angles. For electron 
distributions which are beamed along the field lines, the radiation in the range 9 = 0° 


12 



to 90° is dominated by large angle emission (the flat portion of the cross section) from 
electrons with small pitch angles. As the width of the electron beam decreases the 
variation from 0° to 90° decreases, i.e., contrary to naieve expectation the more highly 
beamed models give more isotropic photon distributions in the observable hemisphere. 
For the same reason, in the isotropic and pancake models the presence of a large 
proportion of electrons at high pitch angles results in a larger degree of anisotropy. For 
emission at lower energies (A: < 100 keV) the directivity and spectral index are nearly 
constant because the cross section does not have the sharp maximum, and the electron 
distributions are broadened due to collisional diffusion. The anisotropy that does exist 
(Petrosian and Leach 1983 and Leach 1984) decreases with increasing a* especially when 
pancake type distributions are included. 

4- Magnetic Field Convergence 

Finally, let us consider the effects of non-uniformities in the magnetic field, 
in particular non-uniformities in the form of loops with field lines which converge 
monotonically with depth s. As mentioned above we consider two different magnetic 
field structures. For models with constant dlnB/ds most of the effects of the 
convergence occur in the corona, while most of the radiation (especially at high energies) 
comes from the footpoints, where because of the small density scale height the effect of 
the convergence is small. The photon and electron distributions are therefore similar 
to those for the uniform field models except that there are relatively more electrons 
“backward” hemisphere with pitch angles in the greater than 7r/2 (i.e., moving to lower 
values of a). This is especially true at high energies where the field convergence is much 
more effective than collisional diffusion in reflecting electrons into the “backwards” 
hemisphere. 

These effects are shown in Figures 7a and 7b where we compare the directivities 
and spectral indices for the standard model (solid lines) and one with isotropic injection 
and constant dlnB/ds with mirror ratio b T = B(st r )/B 0 = 10 (dashed lines). The 
anisotropy and the variation with 6 of the spectral indices decrease with increasing 
convergence. This is a direct result of isotropization of the electron distributions by 
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the magnetic mirroring process. Note that as in the uniform field models, most of the 
angular variation takes place between 60° and 90°. This is a characteristic of all the 
“precipitation” models, i.e., those with field convergence in the corona. 

For models with B(s ) oc (n(s))*, there is no convergence in the corona and 
convergence takes place where the electrons emit most of the radiation. The o’s in 
Figure 7 represent the directivities and spectral indices for a model with mirror ratio 
b T = B(s ma . x )/B 0 = 5. For this type of model the trapping and isotropization become 
more effective at higher energies and we find that unlike the above models there is 
significant difference for all angles, and in some cases (for extreme relativistic energies) 
a significant part of the variation takes place between 0° and 60°. 

5. Combined Effects of a 2 and dhxB/ds 

We have shown the effects of varying a 2 for uniform fields and the effects of field 
convergence for isotropic emission. When both at 2 and B(s) are changed, the effects of 
each variation are combined. It would take a considerable amount of space to present 
the combined effects graphically. Instead, we summarize the results for the constant 
dlnB/ds models and the B oc n* in Tables 1 and 2, repectively. In each table values 
of the total anisotropy from 0° to 90° at 1 MeV, A(fc,0“), and AI\ = r(0°) — T(90 <> ) for 
in the spectral indices T c , r m and I\ are given for 5 different pitch angle distributions. 
Note that this does not show shape of the variation with angle, but the values at 
intermediate angles may be estimated using Figures 6 and 7 as guides. 

In Table 1, the largest variation in directivity and spectral index occurs at larger 
values of field convergence (or b T ) for beamed versus pancake models; it occurs for 
uniform field for highly pancaked models. 

For the B oc n* models, shown in Table 2, the behavior is more complex. In general, 
for a given value of the field convergence (b T ) the disk to limb variation increases from 
left to right (beam to pancake distributions) with the change becoming negligible at 
larger values of b T . For a given value of a 2 the change in variation with b T is largest and 
in the oppisite direction for highly beamed and pancaked models, except for Ar m , which 
has essentially the same variation for all a 2 . 
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For comparison, observations by the GRS on SMM (Vestrand et al. 1987, Bai 1988) 
show (Ar m ) « 0.8 and (log(A(l MeV, 0°))) « —1.2. Inspection of Tables 1 and 2 shows 
that some values of a* and b T can fit this data rather well. Interestingly these include 
both beam and pancake models but none with isotropic injection. In the next paper we 
will make a more detailed comparison with these statistical results and will show that 
beamed models with moderate field convergence are favored. 

6. Breaks in Spectra 

We have seen that the spectral index changes with energy, and the spectral indices 
fit for different energies have different angular dependences. For comparison with the 
data just mentioned, we need to determine the characterisitcs of such spectral breaks. 
Figure 8 shows the difference between the intermediate and low energy spectral indices, 
7m— 7*, for the models with constant dlnB/ds and b T = 1 and 10. A positive (negative) 
value describes a steepening (flattening) of the spectrum. The converging models show 
less steepening for all three injected spectra. For both the uniform and converging 
fields, the change in spectral index is largest for the beam injection, it decreases with 
increasing width ot\ and is smallest for pancake injection. The size of the spectral break 
also depends on the angle, decreasing with 9 for all the models except for the uniform 
field beam injection. For the uniform field case, the spectrum steepens for the beam and 
isotropic models over the whole range of 6. The pancake injection shows some flattening 
at 90°. The size of the break decreases for increasing convergence, and for b T = 25 and 
pancake injection we see flattening for all values of 9. 

The SMM observations mentioned above show that the difference between average 
spectral indices (7*) — (7 m ) is about one with no discemable angular variation. 

None of the models described above can fit this which may indicate that some of our 
assumptions, in particular, a constant power law spectral index for accelerated electrons, 
may be incorrect. We shall discuss this this in detail in Paper III. 

IV. Comparison with Stereoscopic Observations 

We now turn to comparison with the stereoscopic observations of flares made 
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by Kane et al. (1988) leaving a detailed comparison with SMM results to Paper III. 
Kane et al. give photon fluxes and spectral indices for nine large flaxes observed by 
detectors on two satellites, the PVO (for Pioneer Venus Orbiter) and ISEE-3 (for Third 
International Sim Earth Explorer). Ideally one would like use such data to find the 
range of model parameters which fit each observation. This is not possible for the data 
in question because of the large error bars. Instead we compare all of the observations 
with each of the models discussed above (see Tables 2 and 3), to see which of the events 
fit each model. Examples of such cmoparisons are shown in Figures 9 and 10. The 
750 keV fluxes for the two instruments versus angle for the nine events are shown in 
Figure 9 compared with the flux from the model with the beam injection model with 
al — 0.40 and a uniform field. [This choice is motivated by the good agreement that 
this model shows with SMM results to be discussed in Paper III.] The observed PVO 
flux is normalized to the model curve at the appropriate angle, and the ISEE-3 flux is 
shown at its appropriate angle with the same normalization. The error bars are those 
for the flux ratio J(ISEE — 3)/7(PVO) given by Kane et al. The two observed fluxes 
are connected by dashed lines. Five of the flares fit this model within the uncertainties, 
two more fit beamed models with smaller values of a*, and two are not close, especially 
the flare seen near 90° by both instruments. These are the events for which the photon 
flux index can be larger for the larger angle (including uncertainties). In Figure 10, we 
make the same comparison using the spectral indices measured by each instrument. The 
results are similar, but the uncertainties here are larger (to reduce confusion we show 
some average error bars which are two to three times smaller than the actual observed 
error bars), and the observations in some cases do not agree with any model or with 
each other even when the veiwing angles are similar. 

Two of the flares cannot be reconciled with the models, all of which predict a larger 
photon flux mid smaller spectral index with increasing 9. In order to fit such flares we 
must relax our assumptions, in particular the assumption of vertical magnetic fields and 
a low density corona. The geometric factors that have been ignored may be important. 

If the magnetic field is not vertical at the footpoints then the viewing angle and the 
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angle of emission with respect to the field do not correspond to each other and there is 
no reason to expect the photon flux to increase with increasing observation angle. 

V. SUMMARY 

We have described the angular variation of bremsstrahlung radiation with 
observation angle and discussed the effects of different injected distributions and 
magnetic field structures. We find that our models fit some stereoscopic observations of 
real flares, but find it difficult to constrain the flare parameters using these comparisons. 
Taken at face value, this comparison for some flares it seems to imply that our 
assumption that the magnetic field is vertical at the footpoints, and that the electrons 
are accelerated in the corona may not be correct. However, due to the observational 
uncertainties in the data such a conclusion is not warranted yet. In a subsequent paper 
we shall compare with some observations which fit the models much better, and give us 
a reason to believe that our assumptions are correct. 
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Table 1 

Anisotropy D(k ) = A(k, 0°) at 1 MeV, the difference in spectral indices from 0° 
to 90°, Ar = r(0°) — r(90°) for T z , r m and r 7 , for the constant d\nB/ds 
models, for four values of the mirror ratio b T and five values of Beam 
injection is denoted by a (b), and pancake injection is denoted by a (p). 


log [-0(1 

KWo 

MeV)] 

0.04(b) 

0.40(b) 

OO 

0.40(p) 

0.04(p) 

1.0 

-0.93 

-1.29 

-1.72 

-1.91 

-2.14 

2.5 

-1.28 

-1.73 

-1.74 

-1.73 

-1.65 

10.0 

-1.65 

-1.60 

-1.50 

-1.39 

-1.26 

25.0 

-1.57 

-1.39 

-1.27 

-1.08 

-0.95 

Ar, (16 
K Wo 

to 210 keV) 
0.04(b) 

0.40(b) 

oo 

0.40(p) 

0.04(p) 

1.0 

0.39 

0.50 

0.64 

0.73 

0.81 

2.5 

0.48 

0.63 

0.68 

0.70 

0.71 

10.0 

0.61 

0.63 

0.60 

0.56 

0.52 

25.0 

0.61 

0.54 

0.48 

0.38 

0.32 

Ar m (3( 
KWo 

30 keV to 1 MeV) 

0.04(b) 0.40(b) 

oo 

0.40(p) 

0.04(p) 

1.0 

0.26 

0.65 

1.16 

1.33 

1.62 

2.5 

0.70 

1.22 

1.22 

1.21 

1.15 

10.0 

1.15 

1.15 

1.11 

1.06 

0.97 

25.0 

1.11 

1.05 

1.01 

0.92 

0.84 

Ar 7 (2j 
KWo 

5 to 21.2 MeV) 
0.04(b) 

0.40(b) 

00 

0.40(p) 

0.04(p) 

1.0 

0.60 

0.80 

1.47 

1.49 

1.73 

2.5 

1.76 

1.18 

1.02 

0.94 

0.86 

10.0 

1.12 

0.94 

0.88 

0.84 

0.81 

25.0 

0.88 

0.75 

0.71 

0.70 

0.68 
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Table 2 


Same as Table 1 for the models with B(s) oc [n(s)] e . Here there are three values 
of the mirror ratio b T . 


log(£>(l MeV)] 



0.04(b) 

0.40(b) 

OO 

0.40(p) 

0.04(p) 

1.0 

-0.93 

-1.29 

-1.72 

-1.91 

-2.14 

2.5 

-1.46 

-1.90 

-1.97 

-2.00 

-2.00 

5.0 

-1.75 

-1.87 

-1.88 

-1.89 

-1.88 

Ar x (16 

to 210 keV) 





K\«l 

0.04(b) 

0.40(b) 

oo 

0.40(p) 

0.04(p) 

1.0 

0.39 

0.50 

0.64 

0.73 

0.81 

2.5 

0.43 

0.60 

0.66 

0.71 

0.73 

5.0 

0.50 

0.60 

0.63 

0.67 

0.66 

Ar m (31 

)0 keV to 1 MeV) 




KWo 

0.04(b) 

0.40(b) 

00 

0.40(p) 

0.04(p) 

1.0 

0.26 

0.65 

1.16 

1.33 

1.62 

2.5 

0.99 

1.44 

1.50 

1.52 

1.52 

5.0 

1.35 

1.45 

1.45 

1.45 

1.46 

AI\ (2J 

5 to 21.2 MeV) 




b T Wo 

0.04(b) 

0.40(b) 

oo 

0.40(p) 

0.04(p) 

1.0 

0.60 

0.80 

1.47 

1.49 

1.73 

2.5 

2.22 

2.06 

2.02 

2.02 

2.03 

5.0 

2.07 

2.02 

2.02 

2.01 

2.01 
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Figures 


Figure 1: The logarithmic variation with height above the transition region of the 
magnteic field and plasma; nio is the tootal density in units of 10 10 cm -3 , iV 19 is 
the column depth in units of 10 19 cm -2 , X is the temperature, B p and B t are the 
magnetic field for two models: for B p , dlnB/ds = constant, and for B t , B oc n° 2 . 

Figure 2: Fractional emission, (a) r c I(fc,r c ) versus normalized column depth r c , (b) 

I(Jc, h ) versus height h for 5 keV, 16 keV, 300 keV, 1 MeV, and 10.6 MeV emission 
for isotropic injection and uniform field. Each curve is labeledby the photon energy 
in MeV. 

Figure 3: Directivity A(&, 9, r c ) versus observation angle 9 at different depths at 1 MeV 
for the isotropic injection, uniform field model. Each curve is labeled with the 
appropriate value of T c /r)(k). The curve labeled Total is the directivity for the 
whole loop, A (k,9) 

Figure 4: Photon fluxes (multiplied by k 6-1 ) versus energy k. The solid curves labeled 
corona, photosphere and 90 represent J(k,0 , r c ), the emission integrated over 
t c for the corona, photosphere and the whole loop for a limb flare (9 = 90°). 

The dashed lines represent J(k,9), the emission from the whole loop, for 0 = 0°, 
60°, and the curve labeled Total represents J{k ), the emission integrated over all 
angles. 

Figure 5 : Directivity A(Jfc, 9) versus observation angle 9 for the isotropic injection and 
uniform field model. Each curve is labeled by the value of the photon energy in 
MeV. The dahsed line represents the 1 Mev directivity for a model with constant 
density and fully ionized plasma and the o’s represent a model with density 
increasing in depth, but also fully ionized. 

Figure 6: (a) Directivity versus observation angle for emission at 16 keV, 1 MeV and 
10.6 MeV for uniform magnetic field models. The solid lines denote isotropic 
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injection, dashed lines denote beam injection with = 0.04, and the lines 
denote pancake injection with a 2 0 = 0.04. (b) Same as (a) for the spectral indices 
T* (shifted down by 1.0), F m (not shifted) and (shifted up by 1.0). The curves 
are labeled by the range in energy for the spectral indices. 

Figure 7: Same as Figure 6 for isotropic injection models with uniform field (solid lines) 
constant dhxB/ds with b T = 10 (dashed lines) and B <x n 0 2 (o’s) 

Figure 8: The break in the spectral index r m — T, for models with constant dlnB/ds 
and 6=1 and 10. The solid lines denote isotropic injection, dashed lines denote 
beam injection with a* = 0.40, and the lines denote pancake injection with 
a* = 0.40. The upper curve for each case is the uniform field model, the lower one 
the converging field model. 

Figure 9: A comparison of stereoscopic observations of photon flux at 750 keV from the 
PVO and ISEE-3 spacecraft with numerical results. The curve is the directivity 
versus 6 for the model with beam injection of width ol\ = 0.40 and a uniform 
magnetic field. The open circles are the PVO fluxes normalized to the numerical 
results, and the circles with error bars are the normalized values of the ISEE-3 
photon fluxes the two observations for each event are connected by dashed lines. 
The filled circles represent the events for which no models fit. 

Figure 10 : Same as Figure 9 for the spectral indices 7 , and 7 ,. 
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